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Abstract. 

By using most of the present Cosmic Microwave Background (CMB) and Large 
Scale Structure (LSS) measurements and the Big Bang Nucleosynthesis (BBN) 
constraints on the primordial helium abundance, Yp, we set bounds on the radiation 
content of the Universe and neutrino properties. We consider lepton asymmetric 
cosmological models parametrized by the neutrino degeneracy parameter S.^ and the 
variation of the relativistic degrees of freedom, AA^^//, due to possible other physical 
processes that occurred between BBN and structure formation epoch. 
We found that present CMB and LSS data constraints the neutrino degeneracy 
parameter at ^i, < 0.722, implying a lepton asymmetry of the neutrino background 

< 0.614 (2-ct). We also found AN^^f = 0.572t};^^o , the contribution to the 
effective number of relativistic neutrino species N"^^^ — S.OSS^J^yg and a primordial 
helium abundance Yp = 0.249^o Q}g (2-ct errors). 

These results bring an important improvement over the similar ones obtained by using 
WMAP 1-year and older LSS data or the WMAP 3-year data alone and the standard 
primordial helium abundance value Yp = 0.24, relaxing the stringent BBN constraint 
on the neutrino degeneracy parameter (^^ < 0.07). 

We forecast that the CMB temperature and polarization maps observed with high 
angular resolutions and sensitivity by the future Planck Mission will constraint the 
primordial primordial helium abundance at Yp = 0.247±0.002 (2-cr errors) in agreement 
with the most stringent limits on Yp given by the BBN and the neutrino degeneracy 
parameter at < 0.280 (2-it), not excluding the possibility of larger lepton asymmetry. 

This work has been done on behalf of Planck-LFI activities. 



PACS numbers: CMBR theory, dark matter, cosmological neutrinos, gravitational 
lensing 
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1. Introduction 

The radiation budget of the Universe rehes on a strong theoretical prejudice: apart from 
the Cosmic Microwave Background (CMB) photons, the relativistic background would 
consist of neutrinos and of possible contributions from other relativistic relicts. The main 
constraints on the radiation energy density come either from the very early Universe, 
where the radiation was the dominant source of energy, or from the observation of 
cosmological perturbations which carry the information about the time equality between 
matter and radiation. 

In particular, the primordial light element abundance predictions in the standard 
theory of the Big Bang Nucleosynthesis (BBN) [1, 2, 3, 4] depend on the baryon-to- 
photon ratio, tib, and the radiation energy density at the BBN epoch (energy density 
of order MeV"^), usually parametrized by the effective number of relativistic neutrino 
species, N'^^^f. 

Meanwhile, the number of active neutrino flavors have been fixed by boson decay 
width to Ni, = 2.944 ± 0.012 [4] and the combined study of the incomplete neutrino 
decoupling and the QED corrections indicate that the number of relativistic neutrino 
species is N^-^f — 3.046 [5]. Any departure of N^^^ from this last value would be due to 
non-standard neutrino features or to the contribution of other relativistic relics. 
The solar and atmospheric neutrino oscillation experiments [6, 7] indicate the existence 
of non-zero neutrino masses in eV range. 

There are also indications of neutrino oscillations with larger mass-squared 
difference, coming from the short base-line oscillation experiments [8, 9], that can be 
explained by adding one or two sterile neutrinos with eV-scale mass to the standard 
scheme with three active neutrino flavors (see Ref.[10] for a recent analysis). Such 
results have impact on cosmology because sterile neutrinos can contribute to the number 
of relativistic degrees of freedom at the Big Bang Nucleosynthesis [11]. These models 
are subject to strong bounds on the sum of active neutrino masses from the combination 
of various cosmological data sets [12, 13], ruhng out a thermalized sterile neutrino 
component with eV mass [14, 15]. 

However, there is the possibility to accommodate the cosmological observations with 
data from short base-line neutrino oscillation experiments by postulating the existence of 
a sterile neutrino with the mass of few kcV having a phase-space distribution significantly 
suppressed relative to the thermal distribution. 

For both, non-resonant zero lepton number production and enhanced resonant 
production with initial cosmological lepton number, keV sterile neutrinos are produced 
via small mixing angle oscillation conversion of thermal active neutrinos [16]. 
Sterile neutrino with mass of few keV provides also a valuable Dark Matter (DM) 
candidate [17, 18, 19, 20], alleviating the accumulating contradiction between the ACDM 
model predictions on small scales and observations, by smearing out the small scale 
structure. 
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On the other hand, the possible existence of new particles such as axions and gravitons, 
the time variation of the physical constants and other non-standard scenarios (see e.g. 
[21] and references therein) could contribute to the radiation energy density at BBN 
epoch. 

At the same time, more phenomenological extensions to the standard neutrino 
sector have been studied, the most natural being consideration of the leptonic 
asymmetry [22, 23, 24], parametrized by the neutrino degeneracy parameter — /lu/T^o 
[fXu is the neutrino chemical potential and Ti.^ is the present temperature of the neutrino 



Although the standard model predicts the leptonic asymmetry of the same order 
as the baryonic asymmetry, B ~ 10~^°, there are many particle physics scenario 
in which a leptonic asymmetry much larger can be generated [25, 26]. One of the 
cosmological implications of a larger leptonic asymmetry is the possibility to generate 
small baryonic asymmetry of the Universe through the non-perturbative (sphaleron) 
processes [27, 28, 29]. Therefore, distinguishing between a vanishing and non- vanishing 
at the BBN epoch is a crucial test of the standard assumption that sphaleron effects 
equilibrate the cosmic lepton and baryon asymmetries. 

The measured neutrino mixing parameters implies that neutrinos reach the chemical 
equilibrium before BBN [30, 31, 32] so that all neutrino flavors are characterized by the 

same degeneracy parameter, C,u, at this epoch. 

The most important impact of the leptonic asymmetry on BBN is the shift of the beta 
equilibrium between protons and neutrons and the increase of the radiation energy 
density parametrized by: 



The BBN constraints on N^-^-'^ have been recently reanalyzed by comparing the 
theoretical predictions and experimental data on the primordial abundances of light 
elements, using the baryon abundance derived from the WMAP 3-year (WMAP3) CMB 
temperature and polarization measurements [33, 34, 35]: rjB = 6.14 x 10~^°(1.00±0.04). 
In particular, the ^He abundance, Yp, is quite sensitive to the value of N^^^ . 
The analysis of Ref. [36], the conservative error analysis of helium abundance, Yp — 
0.249 ±0.009 [37], yielded to N^^f = ^-^-li (S-fj) in good agreement with the standard 
value, but still leaving some room for non-standard values, while more stringent error 



bars of helium abundance , Yp = 0.2516 ± 0.0011 [38], leaded to A^^// = 3.32l°|| {2-a) 
[39]. 



The stronger constraints on the degeneracy parameter obtained from BBN [40] 
gives —0.04 < ^ < 0.07 (l-o-), adopting the conservative error analysis of Yp of Ref. [37] 
and $, = 0.024 ± 0.0092 (l-a), adopting the more stringent error bars of Yp of Ref. [41]. 

The CMB anisotropics and LSS matter density fluctuations power spectra carry the 
signature of the energy density of the Universe at the time of matter-radiation equality 
(energy density of order eV^), making possible the measurement of A^'^-^'^ through its 



background, T^jT^^h 



(4/ll)V3]. 




2 



4 



(1) 



Constraints on the lepton asymmetry and radiation energy density: Implications for PlanckA 
effects on tlie growtli of cosmological perturbations. 

More effective number of relativistic neutrino species enliances tlie integrated Saclis- 
Woffe effect on tlie CMB power spectrum, leading to a higher first acoustic Doppler 
peak amplitude. Also, the delay of the epoch of matter-radiation equality shifts the LSS 
matter power spectrum turnover position toward larger angular scales, suppressing the 
power at small scales. In particular, for the leptonic asymmetric models, the neutrino 
mass is lighter than in the symmetric case. This leads to changes in neutrino free- 
streaming length and neutrino Jeans mass due to the increase of the neutrino velocity 
dispersion [42, 43]. 

After WMAP3 data release, there arc many works aiming to constrain N^^f from 
cosmological observations [14, 33, 36, 44, 45, 46]. Their results suggest large values for 
A^^'^'^ within 2-0" interval, some of them not including the standard value N^^^ = 3.046 
[14, 33, 36]. Recently Ref. [47] argues that the discrepancies are due to the treatment of 
the scale-dependent biasing in the galaxy power spectrum inferred from the main galaxy 
sample of the Sloan Digital Sky Survey data release 2 (SDSS-DR2) [48, 49] and the large 
fluctuation amplitude reconstructed from the Lyman— a forest data [50] relative to that 
inferred from WMAP3. 

Discrepancies between BBN and cosmological data results on N^^^ was interpreted as 

2-0" evidence of the fact that further relativistic species are produced by particles decay 
between BBN and structure formation [45, 46]. Other theoretical scenarios include the 
violation of the spin-statistics in the neutrino sector [51], the possibility of an extra 
interaction between the dark energy and radiation or dark matter, the existence of a 
Brans-Dicke field which could mimic the effect of adding extra relativistic energy density 
between BBN and structure formation epochs [52] . 

The extra energy density can be split in two distinct uncorrelated contributions, 
first due to net lepton asymmetry of the neutrino background and second due to the 
extra contributions from other unknown processes: 

The aim of this paper is to obtain bounds on the neutrino lepton asymmetry and on the 
extra radiation energy density by using most of the existing CMB and LSS measure- 
ments and self-consistent BBN priors on Yp. We also to compute the sensitivity of the 
future Planck experiment [53] for these parameters testing the restrictions on cosmo- 
logical models with extra relativistic degrees of freedom expected from high precision 
CMB temperature and polarization anisotropy measurements. 

The paper is organized as follows: Section 2 contains a review on the lepton 
asymmetric cosmological models and the BBN theory. Section 3 is devoted to a summary 
data analysis method while in Section 4 we discuss our results. We draw our main 
conclusions in Section 5 
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Figure 1. The CMB temperature (top panels) and polarization (bottom panels) 
percentage differences corresponding to different helium fractions variations AYp with 
respect to the standard value Yp = 0.248 for: ^ 7^ and AA^^°// = (left panels) and 
^^oth ^ = (right panels). All other parameters are fixed to the values of 
our fiducial model. 



2. Leptonic asymmetric cosmological models and the BBN theory 

The density perturbations in leptonic asymmetric cosmological models have been 
discussed in literature [42, 43, 54, 55]. We applied them to modify the Boltzmann Code 

for Anisotropics in the Microwave Background (CAMB) [56, 57, 58] to compute the 
CMB temperature and polarization anisotropics power spectra and LSS matter density 
fluctuations power spectra for the case of three degenerate neutrinos/antineutrinos with 
mass and degeneracy parameter As neutrinos reach their approximate chemical 
potential equilibrium before BBN epoch [30, 31, 32], we consider in our computation 
that all three flavors of neutrinos/antineutrinos have the same degeneracy parameter 
For simplicity, we also consider all there neutrino/antineutrino flavors with the same 
mass m^. 

When the Universe was hot enough, neutrinos and antineutrinos of each flavor behave 
like relativistic particles with Fermi- Dirac phase space distributions: 

1 1 

where E,y = \/ q'^ + a'^m,y is one flavor neutrino/antineutrino energy and q = ap is the 
comoving momentum. Hereafter, a is the cosmological scale factor (oq = 1 today). The 
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mean energy density and pressure of one flavor of massive degenerated neutrinos and 
antineutrinos can be written as: 

P. + Pp = (ksnr q^E^ifM + fM) , (4) 

3(P. + Pp) = {ksT^f ^ ^UM) + fM) ■ (5) 

We modify in CAMB the expressions for the energy density and the pressure in 
the relativistic and non-relativistic hmits for the degenerate case [43] and follow the 
standard procedure to compute the perturbed quantities by expanding the phase space 
distribution function of neutrinos and antineutrinos into homogeneous and perturbed 
inhomogeneous components [58, 59, 60]. Since the gravitational source term in the 
Boltzmann equation is proportional to the logarithmic derivative of the neutrino 
distribution function with respect to comoving momentum, dln(/,^ + /p)/dlng, we also 
modify this term to account for ^ ^ [43, 54]. 

As mentioned in the first section, the BEN theory gives strong constraints on A^^-^-^ 
at this epoch by comparing the measured light element abundance with the theoretical 
predictions. The only free parameter is the baryon to photon ratio, r]B = nb/n^, that is 
obtained from the CMB observation of Q^/i^. 

In particular, the ^He mass fraction,!^, affects the CMB angular power spectra through 
its impact on different evolution phases of the ionization/recombination history [61]. 

We modify the recombination routine recfast [62] of the CAMB code to explicitly 
account for the dependence of Yp on Qbh"^ and on ^N'^^^ as defined in equation (2), as 
previously suggested in Ref . [63] , by adopting the fitting formula [64] : 

/ 8 8 8 8 

\n=l n=l n=l n=l / 

X exp e„x" , (6) 



where x — log]^o(10^°)^ ^^'^ lO^^rj — 273.49^6/1^. The coefficients a„, c„, dn and 
e„ are given in Ref. [64]. The standard prediction of BBN Yp — 0.248 is obtained for 
AA^^-'^-'^ = 0. According to Ref. [64], the accuracy of this fitting formula is better than 
0.05% for the range 5.48 x 10-^° < Vb < 7.12 x 10"^° (0.02 < Q^/i^ ^ o.026) which 
corresponds to the 3-cr interval obtained by WMAP3 and — 3 < AN^^^ < 3. 
Figure 1 presents the CMB temperature and polarization percentage differences 
corresponding to the different variations of the helium fraction, AYp, with respect to 
the standard value Yp = 0.248 obtained for ^ 7^ and AN^^ = and AN^^ 7^ and 
1^ = 0. The impact of the percent change in Yp on the ionization/recombination history 
has a net impact on the CMB temperature and polarization power spectra at percent 
level. 
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3. Analysis 

We use the CoSMoMC Monte Carlo Markov Chain (MCMC) pubhc package [65] 
modified for our extended 6+3 parameter space to sample from the posterior distribution 
giving the following experimental datasets. 

The Cosmic Microwave Background (CMB): We use the WMAP3 data [33, 34, 35] 
complemented with the CMB data from Boomerang [66, 67], ACBAR [68] and CBI [69] 
experiments. 

Large Scale Structure (LSS): The power spectrum of the matter density fluctuations has 
been inferred from the galaxy clustering data of the Sloan digital Sky Survey (SDSS) 
[48, 49, 70, 71] and Two-dcgrcc Field Galaxy Rcdshift Survey (2dFGRS) [72]. 
In particular, the luminous red galaxies (LRG) sample from the SDSS data release 5 
(SDSS-DR5) has more statistical signiflcance [70, 71] than the spectrum retrieved from 
the SDSS main galaxy sample from data release 2 (SDSS-DR2) [48, 49] ehminating 
the existing tension between the power spectra from SDSS-DR2 and 2dFGRS. For 
this reason we consider in our analysis the matter power spectra from SDSS-LRG and 
2dFGRS. We consider SDSS-LRG data up to kmax ^ 0.2h Mpc'^ and the 2dFGRS data 
up to kmax — 0.14h Mpc^^. We apply the corrections due to the non-linearity behavior 
and scale dependent bias as indicated in Rcf.[70], connecting the linear matter power 
spectrum, Piin{k), and the galaxy power spectrum, Pgai{k), by: 

P^aiik) = b^^-^^^PUk) , (7) 
where the free parameters b and Qni are marginalized. 

Type la Supernovae (SNIa): We also use the luminosity distance measurements of 
distant Type la supernovae obtained by Supernova Legacy Survey (SNLS) [73] and the 
Hubble Space Telescope [74] . 

Hubble Space Telescope key project (HST): We impose priors on the Hubble constant 

i^o = 72 ± 8 km s^-'^Mpc^-'^ from HST key project [75]. 

BEN constraints on Yp-. We use the BEN constraints on Yp as given in equation (7), 
allowing Qbh^ and AN^^^ to span the following ranges: 0.02 < Vti^h"^ < 0.026 and 
-3 < AN^ff < 3. 

Hereafter, we wiU denote WMAP3+SDSS-DR5+2dFGRS+SNIa+HST+BBN data set 
as WMAP3+AU. 

We perform our analysis in the framework of the extended ACDM cosmological 
model described by 6 -|- 3 free parameters: 

standard 

Here Qbh"^ and Qcdmh"^ are the baryonic and cold dark matter energy density parameters, 
6s is the ratio of the sound horizon distance to the angular diameter distance, r 
is the reionization optical depth, Ug is the scalar spectral index of the primordial 
density perturbation power spectrum and Ag is its amplitude at the pivot scale 
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Table 1. The free parameters of our model, their fiducial values used to generate the 
PLANCK-like simulated power spectra and the prior ranges adopted in the analysis. 



Parameter Fiducial value Prior range 

YlJ? 0.022 0.005 ^ 0.1 

Qcdmh^ 0.105 0.01 ^ 0.5 

Os 1.04 0.5 ^ 5 

r 0.09 0.01 ^ 0.3 

Us 0.95 0.5 ^ 1.3 

ln[10i%] 3 2.7 ^ 4 

U 0.05 ^ 0.5 

0^4 

AiV;//^ 0.046 -3 ^3 

Yp 0.248 



— 0.05 hMpc~^. The additional three parameters denote the neutrino energy density 
fraction fi^, the neutrino degeneracy parameter and the extra contributions from 
other unknown processes AN^^fl . Table 1 presents the parameters of our model, their 
fiducial values used to generate the PLANCK-like simulated power spectra and the prior 
ranges adopted in the analysis. 

For the forecast from PLANCK-like simulated data we use the CMB temperature (T) 
and polarization (P) power spectra of our fiducial cosmological model and the expected 
experimental characteristics of the Planck frequency channels presented in Table 2 
[53]. For each frequency channel we consider an homogeneous detector noise with the 
power spectrum: 

where u is the frequency of the channel, 6b is the FWHM of the beam and Ac are the 
corresponding sensitivities per pixel. The global noise of the experiment is obtained as: 



(10) 



In order to interpret the likelihood function, C{Q), as probability density we as- 
sume uniform prior probability on the parameters (i.e. will assume that all values 
of parameters are equally probable). For each parameter we compute the cumulative 
distribution function C(6) = C(Q)dQ/ /if"""' £(i0 and quote as upper and lower 
intervals at 2-(j the values at which C{6) is 0.95 and 0.05 respectively. For the case 
when C{Q) is zero and has a positive values (i.e. f,^, the absolute value of we 
quote only the upper limit at 2-0". 
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Table 2. The expected experimental characteristics for the Planck frequency 
channels considered in the paper. At and Ap are the sensitivities per pixel for 
temperature and polarization maps. 
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Figure 2. The marginalized posterior likelihood probabilities of the main cosmological 
parameters obtained for: WMAP3+A11 with AN^^j^ = prior (black lines) and 
WMAP3+AU (red lines) and Planck (green lines) without priors on AN^H . 

4. Results 

We start by making a consistency check, verifying that by using WMAP3+A11 data and 
imposing = 0, AN^ff^ = and Yp = 0.248 priors we obtain results in agreement with 
the ones obtained by WMAP collaboration (Tables 5 and 6 from Ref. [33]). 

In order to understand how the extra relativistic energy density and the leptonic 
asymmetry affect the determination of other cosmological parameters, we compute first 
the likelihood functions for WMAP3+A11 by imposing AN^^if = prior. We then 
extend our computation over the whole parameter space for WMAP3+A11 and Planck- 
like simulated data. In Figure 2 we compare the marginalized likelihood probabilities 
obtained for WMAP3+A11 with AA^^f/ = prior with those obtained for WMAP3+A11 
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Figure 3. The 2D marginalized 1-a and 2-a allowed regions in to - \^\ and to - AN^^j^ 
planes. The black curves correspond to WMAP3+A11 with AiV^// ~ prior. The red 
and green curves correspond respectively to WMAP3+A11 and PLANCK-like simulated 
data, without priors on ATV^//. 



and Planck without priors on AN^ijj . The main effect of including the contribution 
of the extra relativistic energy density is the change in the age of the Universe (and in 
the Hubble expansion rate) from to = 13.81 ± 0.26 GYrs to to = l3A2t\i2 GYrs {2-a 
errors), effect that is mostly driven by the increased degeneracy between matter and 
radiation energy densities. 

We present in Figure 3 the 2D marginalized 1-a and 2-a allowed regions in to - \^\ and 
to - '^^oth pla-^iss showing this effect. The 2-a confidence region for the additional 
number of relativistic relicts is -1.207 < AN^Jf[ < 2.572 for WMAP3+A11 and 
-0.226 < AN^th < 0.236 for PLANCK-like simulated data. The negative values of 
^■^oth lowering the amplitude of the CMB and LSS power spectra that can be 
compensated by larger values of the degeneracy parameter. 

Figure 4 presents the 2D marginalized 1-a and 2-a allowed regions in a^ - Ug plane and 
in - r plane obtained for WMAP3+A11 with AA^^^f/ = prior and WMAP3+A11 
and PLANCK-like simulated data without priors on AN^^f^ . It is evident from these 
plots that the inclusion of an additional number of relativistic relicts increases the 
degeneracy between the cosmological parameters controlling the CMB and LSS power 
spectra amplitudes. 

In Figure 5 we compare the marginalized likelihood probabilities of the neutrino 
parameters and the helium mass fraction obtained for WMAP3+A11 with AA'^f/ = 
prior with those obtained for WMAP3+A11 and Planck without priors on AN^j^^ . The 
expectation values and the corresponding errors or the upper limits {2-a) are presented 
in Table 3. 

We show in Figure 6 that the precise measurements of the CMB temperature and 
polarization power spectra from Planck will reduce the degeneracy between and 
AN^fff , allowing better constrains of the scenarios involving additional number of rela- 
tivistic relicts and leptonic asymmetry. 
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Figure 4. The 2D marginalized 1-a and 2-a allowed regions in erg - Ug plane (left 
panel) and the same confidence regions in Us - t plane (right panel) obtained for 
WMAP3+A11 with AiV^^/ = prior (black Unes) and WMAP3+AU (red Hnes) and 
PLANCK-like simulated data (green lines) without priors on AN'^H . 







Figure 5. The marginalized posterior likelihood probabilities of the neutrino 
parameters and helium mass fraction obtained for: WMAP3+A11 with AA^' 



th 







prior (black lines) and WMAP3+A11 (red lines) and Planck (green lines) without 
priors on AN^^jf . 



Constraints on the lepton asymmetry and radiation energy density: Implications for PhANCKl2 

Table 3. Constraints on neutrino parameters and helium mass fraction. The errors 
and the upper limits are at 2-a. 





WMAP3+A11 


WMAP3+A11 


Planck 


Parameter 


^KL' = 




^Kll + 




< 0.033 


< 0.037 


< 0.036 




< 0.590 


< 0.722 


< 0.280 




< 0.833 


< 1.243 


< 0.158 




< 0.474 


< 0.614 


< 0.179 


Neff 


< 3.873 


"-••^ ' ^-1.780 

3.058lJ;?78 


0.008+°-! 
2.920^°:^?^ 




< 0.249 


n 249+0 014 

'-'■^^^-0.016 


0.247 ±0.002 




Figure 6. The 2D marginalized 1-ct and 2-ct allowed regions in |^| - AN^^j^ and 
ICI - '710 planes. The black curves correspond to WMAP3+A11 with IS.N'^^ = prior. 
The red and green curves correspond respectively to WMAP3+A11 and PLANCK-like 
simulated data, without priors on AN^^fj^ . 

5. Conclusions 

In this paper, we set bounds on tlie radiation content of the Universe and neutrino 
properties by using most of the present Cosmic Microwave Background (CMB) and 
Large Scale Structure (LSS) measurements. We also take into account the Big Bang 
Nucleosynthesis (BBN) constraints on the primordial helium abundance, Yp, which prove 
to be important in the estimation of cosmological parameters from future Planck 
data, both in non-degenerate and degenerate BBN models; the importance of the self- 
consistent BBN prior on Yp was also emphasized in two other recent analysis [76, 77]. 
We consider lepton asymmetric cosmological models parametrized by the neutrino 
degeneracy parameter and the variation of the relativistic degrees of freedom, AN^j^^ , 
due to possible other physical processes that occurred between BBN and structure 
formation epoch. 
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We found that present CMB and LSS data together with BBN prior on the primordial 
hehum abundance (1^) constraints the neutrino degeneracy parameter at ^i, < 0.722, 
leading to a lepton asymmetric neutrino background of C^, < 0.614 {2-a). We also found 
^^oth ~ 0-^''^2ltj;'78o ) the contribution to the effective number of re- 
lativistic neutrino species A^^'^-'^ = S.OSSlj^'^yg and a primordial helium abundance 
Yp — 0.249lQ;oJg (2-(j errors). These values represent important improvements over 
the similar results obtained by using WMAP 1-year together with older LSS data [78] 
or the WMAP3 data alone [42] and the standard primordial helium abundance value 
Yp — 0.24, relaxing the stringent BBN constraint on the neutrino degeneracy parameter 
ii. < 0.07). 

We observe that, when using WMAP3+A11 data, an additional number of relativis- 
tic relicts brings a substantial degeneracy in the Qm — plane and a weaker constraint 
on the age of the Universe; the same degeneracy occurs also between other cosmological 
parameters under the same conditions. We therefore conclude that the present cosmo- 
logical data do not favor the variation of the relativistic degrees of freedom, AN^^^ , due 
to other possible physical processes that occurred between BBN and matter-radiation 
equality epoch. 

We forecast that the CMB temperature and polarization maps observed with high 
angular resolutions and sensitivity by the future Planck Mission will constraint the 
primordial primordial helium abundance at 1^ = 0.247± 0.002 (2-a errors) in agreement 
with the most stringent limits on Yp given by the BBN and the neutrino degeneracy 
parameter at < 0.280 (2-cr), allowing larger lepton asymmetry models. Also, they 
will reduce the degeneracy between |^| and AN^^ allowing a better distinction between 
extra radiation energy density coming from an additional number of relativistic relicts 
and from a lepton asymmetric neutrino background. 
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